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Abstract. A new method is presented for recovering the abundances of Damped Ly a systems (DLAs) taking 
into account the effects of dust depletion. For the first time, possible variations of the chemical composition of 
the dust are taken into account in estimating the depletions. No prior assumptions on the extinction properties 
of the dust are required. The method requires a set of abundances measured in the gas and a set of parameters 
describing the chemical properties of the dust. A large subset of these parameters is determined from interstellar 
observations; the others are free parameters for which an educated guess can be made. The method is able to 
qq ■ recover the abundances of the SMC starting from SMC interstellar measurements apparently discrepant from the 

| stellar ones. Application of the method to 22 DLAs with available [Fe/H] and [Si/Fe] measurements gives the 

following results: (1) the mean metallicity of the corrected data is < [Fe/H] > ~ —1.0 dex, about 0.5 dex higher 
\^ ' than that of the original data; (2) the slope of the [Fe/H] versus redshift relation is steeper for the corrected data 

(m ~ —0.3 dex) than for the original ones (m ~ —0.2 dex); (3) the corrected [Si/Fe] ratios are less enhanced, on 
average, than those found in Galactic stars of similar, low metallicity; (4) a decrease of the [Si/Fe] versus [Fe/H] 
ratios, expected by "time delay" models of chemical evolution, is found for the corrected data; (5) the [Si/Fe] 
ratios tend to increase with redshift once corrected; (6) consistency between [Si/Fe] and [S/Zn] measurements, 
O i 1 two independent estimators of the a/Fe ratio, is found only for the corrected abundances. 
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1. Introduction DLA abundances can be severely affected by the effects 

of dust. Evidence of dust in DLAs comes from extinction 

High resolution spectroscopy of QSO absorption line sys- and abundance studies . The re ddening of QSOs lying be- 

; terns is a very powerful tool for probing structures in the hind DLA absorption is a clear sign of dust extinction (Pei 

, early universe, including galaxies in the faint end of the gt ^ lggl; pd & ^ lgg5) Eyidence of elemental dep l e - 

luminosity function. With this technique, elemental abun- Uon similar tQ ^ observed in the nearby ISM (Savage & 

dances can be measured in galaxies up to redshift z « 5 Sembach 1996) was first reported by Pet tini et al. (1994) 

and used to probe the early stages of chemical enrichment. frQm ^ behaviour of thc Zn/Cr ratio in DLAs . The en _ 

In fact, abundance measurements have been obtained for a hancement of the Zn/Cr ratio relative to the solar value 

relatively large number of damped Ly a systems (DLAs), wag interpreted as being due to a d i ffere nt degree of incor- 

the class of QSO absorbers most clearly associated with poration into dugt grains of zinCj a volatUe element; and 

intervening galaxies (Lu et al. 1996, 1998; Pettini et al. chromi a re f ra ctory element. Similar results have been 

1997, 1999, 2000; Prochaska & Wolfe 1999; Centurion et found for othep irg Qf elements ^ „ d of 

al. 1998, 2000; Molaro et al. 2000; Petitjean et al. 2000; affinity ^ dugtj e g Zn/Fe Recent abundance studies 

Prochaska et al. 2001). give furtner support to the presence of dust depletion in 

Abundance measurements can be rather accurate in DLAs ( Hou et al 2 001, Prochaska & Wolfe 2002), indicat- 

DLAs: column densities of H I and of metal species can ing that abundance s taken at face value are of little use in 

easily be determined with errors < 0.1 dex or even < 0.1 understanding the chemical properties of DLA galaxies, 
dex; ionization corrections are generally negligible for 

most of the elements considered (Vladilo et al. 2001; see Various approaches have been used to circumvent 

also Izotov et al. 2001). However, the final accuracy of the problem of dust depletion in studies of DLAs abun- 

dances. One way is to focus on elements essentially undc- 

Send offprint requests to: G. Vladilo pleted, such as N (Lu et al. 1998, Centurion et al. 1998), 
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(Molaro et al. 2000, Dessauges-Zavadsky et al. 2001; 
Levshakov et al. 2002), S (Centurion et al. 2000), and Zn 
(Pettini et al. 1997, 1999; Prochaska & Wolfe 1999, 2001; 
Vladilo 2000). Another is to study DLAs with low dust 
content (Pettini et al. 2000; Molaro et al. 2000, Lopez et 
al. 2002). Although these studies do not depend, essen- 
tially, on the properties of the dust, their application is 
limited to certain elements and to particular sub-samples 
of DLAs. The only way to perform a general study of 
DLAs abundances is to quantify the effects of dust deple- 
tion. 

By comparing the abundance ratios observed in DLAs 
with those measured in the interstellar gas of the Galaxy 
it is possible to derive some indication on the depletions of 
DLAs. Studies of this kind show a similarity to the deple- 
tions typical of Galactic warm gas (Lauroesch et al. 1996; 
Kulkarni et al. 1997; Savaglio et al. 2000). In this type 
of work, two different types of intrinsic DLA abundances 
have been considered: solar abundances and abundances 
typical of metal-poor stars in the Milky Way. Clearly, it 
is rather difficult to disentangle the intrinsic abundance 
patterns of DLAs with this type of approach. 

A method for recovering the intrinsic abundances in 
individual systems was presented by Vladilo (1998, here- 
after Paper I). In that work the chemical composition 
and extinction properties of the dust were assumed to 
be the same as those of the dust in Galactic warm gas. 
The amount of dust extinction for a given level of metal- 
licity, i.e. the dust-to-metals ratio, was allowed to vary 
among DLAs. The dust-to-metals ratio of individual DLAs 
was estimated assuming that zinc tracks the iron-peak el- 
ements and that the observed overabundances of Zn/Cr or 
Zn/Fe are due entirely to differential depletion. A limita- 
tion of that method is that the dust composition may actu- 
ally vary in different environments, even within the Milky 
Way (Savage & Sembach 1996). In addition, some recent 
investigations indicate that zinc may not exactly follow 
other iron-peak elements (Primas et al. 2000; Umeda & 
Nomoto, 2001). According to Prochaska & Wolfe (2002), 
the assumptions adopted in Paper I may force the dust- 
corrected abundances to become closer to solar values. To 
overcome these limitations, a new procedure for estimat- 
ing dust depletion in DLAs is presented in this paper. 
The equations of the procedure are based on a new ex- 
pression for scaling dust depletions up or down accord- 
ing to changes of the physical/chemical conditions of the 
interstellar gas (Vladilo 2002, hereafter Paper II). This 
expression allows the dust chemical composition to vary 
according to changes of the dust-to-metals ratio and of 
the intrinsic abundances of the medium. With this new 
formulation, deviations of the intrinsic Zn/Fe ratios from 
solar values can be consistently accounted for. At variance 
with Paper I, the present method does not make any as- 
sumption concerning the extinction properties of the dust 
in DLAs. The method is presented in Section 2, some ex- 
amples of application are described in Section 3, and the 
results are discussed and summarized in Section 4. 



2. The method 

The starting point of the method is the expression for 
scaling depletions up or down derived in Paper II: 



(1) 



where fx = iVx,dust/A r x,medium is the fraction in dust of 
the element X; the index "medium" refers to the total- 
ity of the atoms (gas and dust); the index i indicates a 
phase of the Galactic ISM for which fx,i can be deter- 
mined empirically; the index j the ISM of an external 
galaxy (for instance, a DLA system); Y is an element 
used as a reference for expressing relative abundances. 
The scaling law ([j]) is derived assuming that the abun- 
dances in the dust, px = (Nx/Ny)dust, are a function 
of the dust-to-metals ratio, r, and of the intrinsic abun- 
dances, ax = (Nx/N Y ) medium, i.e. 



Px =Px(r, a x ) ■ 



(2) 



The dust-to-metals ratio is defined as the fraction in dust 
of the reference element, i.e. r = fy. At variance with 
Paper I, the extinction properties of dust grains do not 
appear in this definition. In Eq. ([l]) the intrinsic abun- 
dances (Nx/Ny) medium are accounted for in the term 
[£] = log(A r x/A r Y)medium - log(X/Y) Q . The parameters 
rjx and ex are dimensionless derivatives of px with respect 
to r and to ax, respectively. Given the function (||), the 
scaling law (|l|) is always correct for infinitesimal changes 
of r and ax- For discrete variations of these parameters, 
the law is valid if rjx and ex are constant. Comparison 
with the observations indicates that a set of constant pa- 
rameters is able to model all the typical depletion patterns 
of the Galactic ISM by only varying r in Eq. (P . By allow- 
ing both r and ax to vary, the few extragalactic patterns 
observed so far[] can be modeled as well. We refer to Paper 
II for more details. 

The depletion patterns that can be successfully mod- 
eled by Eq. ([!]) originate in different types of interstellar 
gas, characterized by a variety of physical conditions (e.g. 
cold gas vs. warm gas), metallicities (e.g. 1 vs. ~ 0.25 
of the solar level), locations (e.g., Galactic halo vs. disk), 
and extinction properties (e.g., strong vs. weak 217.5 nm 
emission bump). For these reasons we feel justified in ap- 
plying the same scaling law to the ISM of distant galaxies, 
including DLAs. 

2.1. Equations of the method 

The relation between the observed and the intrinsic abun- 
dances can be derived from the Eq. ([!]) above and from 
Eq. (4) of Paper I: 
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1 Only some interstellar lines of sight to the Magellanic 
Clouds are currently available. 
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where g = Vj/ri is the dust-to-metals ratio normalized 
to r,;, [X/H] bs the abundance observed in the gas phase, 
pC/H|j the intrinsic abundance and [§] = log(Nx/Nn) — 
log(X/Y) . 

In order to recover the intrinsic abundances we first 
need to determine g. To do this, we use the expression 



Table 1. Adopted interstellar parameters 



0, (4) 



which can be derived by applying Eq. (^|) to the generic 
element X and to the reference element Y. In practice, the 
solution of Eq. requires a measurement of [X/Y] Q b s and 
an educated guess of [X/Y]j. 

Once g is known, we can determine the relative abun- 
dance in the medium [X/Y]j from the expression 



10 



which can be derived by rearranging Eq. (||). Clearly, this 
second expression can be applied to elements X different 
from the one used in Eq. (|4|). 

Eqs. (Q) and (||) can be solved analytically only in par- 
ticular cases, depending on the values of the parameters 
?7x and ex- These equations have the general form 



x ± ax 



b = 



(6) 



and can be solved by iterations given the coefficients a and 
b and the exponent a. In general, the uniqueness of the so- 
lution is not guaranteed. If multiple solutions are present, 
one needs to restrict the search for the roots within a 
limited interval of x. This interval has to be justified on 
physical grounds. 

In the case of Eq. ^ x = g and a = l + r/x- The expo- 
nent a can attain values higher than unity (see adopted 
values of rjx in Table 1). A unique solution is found search- 
ing for the roots in the interval < g < gn m , where 
Qum = 1-0649 is the dust-to-metals ratio corresponding 
to a fraction of iron in dust of 100%. Values outside of 
this interval are obviously unphysical. 

In the case of Eq. (0) x = loM and a — ex- As 
explained below, a choice of the exponent in the range 
< ex < 1 allows us to probe two extreme possibilities 
concerning the dependence of the dust composition on the 
composition of the medium. For these values of a the solu- 
tions of Eq. (||) are unique and there is no need to restrict 
the search for the roots to a particular interval of intrinsic 

abundances x = ioffl. 



2.2. Input parameters 

To apply the procedure various sets of interstellar param- 
eters are required. The fractions in dust fx.i and the pa- 
rameters rjx can be determined from Galactic interstellar 
studies, as shown in Paper II. In Table 1 we list the set 



Element 



7?X 



Si 


0.691 ±0.069 


±4.68 ± 1.77 


Mg 


0.715 ±0.056 


±3.60 ±1.50 


Mn 


0.880 ±0.019 


±0.53 ±0.41 


Cr 


0.920 ±0.010 


±0.37 ±0.20 


Fe 


0.939 ± 0.004 


±0.00 ±0.07 


Ni 


0.941 ± 0.005 


±0.01 ±0.11 


Zn 


0.587 ± 0.048 6 


±3.97 ± 1.56 b 



a The r/x. parameters have been determined by comparing the fractions 
in dust of the warm gas with those of the cold gas in the Galactic disk. 
b Warm disk fraction in dust not available for zinc; the adopted values 
are based on the observed zinc depletions in the cold-idisk and in the 
warm disk+halo gas scaled up or down with relation (111). 



of fx and ijx parameters that allow one to model simul- 
taneously all types of Galactic depletion patterns by only 
varying g = Tj/ri. The Galactic depletion patterns include 
the cold disk, warm disk, warm disk+halo, and warm halo 
ones, as defined by Savage & Sembach (1996)0. 

In principle, the £x parameters can be determined 
from extragalactic interstellar observations. The prelim- 
inary study of the Small Magellanic Cloud described in 
Paper II indicates that the observed depletions can be re- 
produced for ex — 1) or somewhat below unity. However, 
this result could be specific to the few extragalactic lines of 
sight investigated. Therefore, in the present work we prefer 
to consider ex as a free parameter. Luckily, an educated 
guess of ex can be made considering the physical meaning 
of this parameter. As shown in Paper II, ~ £x^^. In 
other words, £x indicates the sensitivity o fthe abundance 
ratio X/Y in the dust to changes of the same ratio in the 
medium. The case £x = 1 is what one may expect on 
intuitive grounds: the more abundant an element is, the 
more atoms of that element will be available for being in- 
corporated into the dust. However, given the complexity 
of the processes involving dust formation, accretion and 
destruction, the composition of the medium might play a 
secondary role in determining the composition of the dust, 
at least for some elements, in which case £x ~ 0. Values 
£x > 1 do not have a special appeal from the intuitive 
point of view and are not supported by the interstellar 
study of the SMC carried out in Paper II. Based on these 
considerations, we adopt here £x = 1 or £x = as two lim- 
iting possibilities. The diversity (similarity) between the 
results obtained in these two cases indicates the weakness 
(robustness) of the method due to the poor knowledge of 
the function px = Px(ax)- Given the form of the equations 
(H) i @ an d (H) ! the uncertainty of £x will significantly af- 



2 The depletions from Savage & Sembach (1996) have been 
updated according to revised values of Mgn (Welty et al. 1999 
and refs. therein) and Nill (Fedchak & Lawler 1999; Howk et 
al. 1999) oscillator strengths. Zinc has been added to the orig- 
inal list by including measurements of the same lines of sight 
selected by Savage & Sembach (Roth & Blades 1995). Solar 
reference values from Anders & Grevesse (1989) are adopted 
consistently in the present paper. 
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Table 2. Sets of [Zn/Fe].,- and ex parameters 



Set 


fZn/Fcl, 






SOO 


+0.0 


o 


o 


S10 


+0.0 


1 





SOI 


+0.0 





1 


Sll 


+0.0 


1 


1 


E00 


+0.1 








E10 


+0.1 


1 





E01 


+0.1 





1 


Ell 


+0.1 


1 


1 



feet the results only if the intrinsic abundances deviate 
significantly from solar ratios. 



2.2.1. The intrinsic [Zn/Fe] ratio 

Interstellar studies indicate that the fraction of zinc in 
dust is relatively small while iron is heavily depleted 
(Savage & Sembach 1996). On the other hand, studies 
of Galactic stars indicate that the zinc abundance nearly 
tracks that of iron-peak elements down to metallicities 
[Fe/H] ~ -2.2 dex (Sneden et al. 1991). Should zinc track 
iron in DLAs precisely (i.e., [Zn/Fe] j =0), the [Zn/Fc] 
ratio would be a perfect indicator of depletion and an 
optimal choice for determining the dust-to-metals ratio 
with Eq. (Q). However, Galactic stars show a significant 
enhancement of the ratio [Zn/Fe] at metallicities [Fe/H] 
< —2.2 dex (Primas et al. 2000). This result suggests a 
difference in the nucleosynthetic origin of zinc and iron 
(Umeda & Nomoto 2002) and makes the use of the [Zn/Fc] 
ratio as an indicator of depletion less reliable (Prochaska 
k Wolfe 2002). The metallicities of DLAs are generally 
higher than —2.2 dex. In this range the [Zn/Fe] ratio in 
Galactic stars is consistent with the solar value within 
the errors, with a scatter of about ±0.1 dex at la level 
and tentative evidence for an enhancement of c± +0.1 
dex. Therefore, we may expect that the enhancement of 
the intrinsic ratio [Zn/Fe]j in DLAs is little or negligible 
in the same range of metallicities. The DLAs with the 
lowest metallicity and dust content should reveal an in- 
trinsic enhancement, if present. The best example at our 
disposal, the DLA at z — 3.39 towards Q0000-26, with 
[Zn/H] ~ -2.1 dex (Molaro et al. 2000), shows very little 
evidence for [Zn/Fe] enhancement. In the rest of this pa- 
per we consider [Zn/Fe],,- as an input parameter to which 
we assign values 0.0 or +0.1 dex in the solution of our 
problem. 

2.3. Abundances of individual velocity components 

The absorption spectrum of DLAs is often characterized 
by the presence of multiple components. These may be 
associated with parcels of gas located in different regions 
of the intervening galaxy. The method presented here 
can be applied to the analysis of individual components. 



However, while the velocity components of metal lines may 
be resolved in high resolution spectra, those of Hi lines 
are not. This implies that Eqs. (ji|) and (0) can be ap- 
plied to individual components, but Eqs. (g) cannot. In 
other words, relative abundances and dust-to-metals ra- 
tios can be obtained for individual components. The ab- 
solute abundances can be corrected for dust effects only as 
far as the integrated absorption is concerned (i.e. the total 
column densities inside the DLA). This contraint is set by 
the nature of the measurements and not by a theoretical 
limitation of the method. 

3. Applications of the method 

Some examples of application of the method are presented 
in this section. First, the method is tested in a nearby 
galaxy, namely the Small Magellanic Cloud. The SMC 
is the only external galaxy with several accurate inter- 
stellar measurements available and with intrinsic abun- 
dances known from stellar observations. We show that the 
method is successful in recovering the intrinsic SMC abun- 
dances, a result that gives us more confidence in applying 
the same procedure in high redshift DLAs with unknown 
stellar abundances. In Section 3.2 the method is used to 
correct [Fe/H] metallicities and [Si/Fe] ratios in DLAs. 

In all the applications presented here Fe is adopted as 
the reference element Y. The Zn/Fe ratio is used to deter- 
mine the normalized dust-to-metals ratio g with Eq. (|J). 
Errors in g are estimated by taking into account the mea- 
surement errors of the Zn n and Fe n column densities. The 
dust-corrected abundances are obtained by solving Eq. 
analytically for ex = and = 1. The errors of the dust- 
corrected abundances shown in Table 3 and in the figures 
are derived from the propagation of the g errors and of the 
[X/Fe] bs or [Fe/H] D b s errors. The adopted parameters fx 
and r/x derived from interstellar studies are listed in Table 
1. The sets of educated guesses of the [Zn/Fe]j and ex pa- 
rameters are shown in Table 2. Since we are interested in 
deriving the absolute abundance of the reference element, 
[Fe/H], the analysis of individual velocity components is 
not considered (see Section 2.5). 

3.1. SMC interstellar lines of sight 

The interstellar lines of sight of the SMC considered here 
are those towards the stars Sk 78, Sk 108, and Sk 155, 
for which accurate abundance measurements are available 
(see Paper II). In Fig. |l]we plot the [Si/Fe] versus [Fe/H] 
results for these three SMC lines of sight. For the sake of 
clarity, only one set of Table 2 input parameters is con- 
sidered in the figure. The uncorrected abundances (filled 
symbols) suggest the existence of a spread of the [Fe/H] 
metallicity and of the [Si/Fe] ratio among the different 
lines of sight. However, none of them is in agreement with 
the abundances representative of the SMC (box in the fig- 
ure). The corrected ratios (open symbols) yield instead 
consistent results, all in good agreement with the intrin- 
sic SMC ratios. Since the three different lines of sight are 
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Fig. 1. Interstellar [Si/Fe] versus [Fe/H] in three lines of 
sight of the SMC (squares: Sk 78; diamonds: Sk 108; cir- 
cles: Sk 155). Filled symbols: observed values taken from 
Welty et al. (2001). Open symbols: values corrected for 
dust effects using the set of parameters Ell. Box: range 
of [Si/Fe] vs. [Fe/H] values obtained from SMC stellar ob- 
servations (Russell & Dopita 1992). 



As a first application of the new method to DLAs we 
present here a list of dust-corrected [Fe/H] and [Si/Fe] 
abundances. Only DLAs with good quality measurements 
of Hi, Sin, Fen, and Znii have been considered. The Zn 
data are required to derive the dust-to-metals ratio g. The 
resulting sample of 22 systems is shown in Table 3 together 
with the relevant references for the adopted abundances. 
The [Fe/H] and [Si/Fe] measurements with and without 
dust correction are compared in the table. Four sets of cor- 
rected abundances are shown, corresponding to the sets of 
input parameters S00, Sll, E00, and Ell in Table 2. The 
results are not strongly affected by the choice of these pa- 
rameters. 

The values of g obtained from the [Zn/Fe] data in- 
dicate that the level of depletion is always smaller than 
that found in the cold gas of the Galactic disk. In fact, 
the median g value of the sample is 0.74 and 0.61 for 
[Zn/Fe]j = 0.0 and +0.1 dex, respectively. These val- 
ues of g imply that the iron fraction in dust in DLAs 
is lower, on average, than that in Galactic warm halo gas, 
Qwh = 0.84 (Paper II). The broad similarity between de- 
pletion in DLAs and depictions in warm halo gas found in 
previous work (Lauroesch et al. 1996; Savaglio et al. 2000) 
is therefore confirmed by the present analysis. In the rest 
of this section we present some general conclusions derived 
from the dust-corrected sample of Table 3. 



treated independently, each with its own value of the g pa- 
rameter, it is remarkable that they give consistent results 
only when the dust correction is applied. 

Similar results are obtained from the study of the 
[Mg/Fe], [S/Fe], [Cr/Fe], [Mn/Fe], and [Ni/Fe] ratios in 
the same SMC lines of sight. These different ratios are 
corrected independently of each other, since each element 
X has its own value of r/x- Therefore, it is remarkable that 
each one of them matches the corresponding SMC [X/Fc] 
ratio only when the dust correction is applied. The results 
are not significantly affected by the choice of the input 
parameters listed in Table 2. 

In summary, the interstellar abundances not corrected 
for dust depletion do not match the intrinsic SMC abun- 
dances. The corrected abundances yield consistent results 
for different elements and different lines of sight, all in 
agreement with the intrinsic SMC abundances. Most of 
the results are modestly affected by the choice of the 
[Zn/Fe]j and £x parameters. 

3.2. Dust-corrected abundances in DLAs 

Iron is often used as a tracer of metallicity in DLAs 
because many Fe II transitions with different oscillator 
strengths can be detected in a broad range of absorption 
redshift. Silicon is the only ev-capture element which has 
been extensively observed in DLAs, with over 30 mea- 
surements available. Therefore, the Si/Fe ratio is the only 
abundance ratio of a-capture over iron-peak elements for 
which a statistically significant sample can be analyised. 



3.2.1. Metallicity 

As in the case of SMC interstellar data, also in DLAs 
the use of [Fe/H] without dust correction can lead to a 
severe underestimate of the real metallicity of the sys- 
tem. As shown in Table 3, the difference between be- 
tween original and corrected measurements can be as 
high as 1 dex in some cases. The statistical properties 
of the sample are severely affected by this effect. For in- 
stance, the mean metallicity of the original measurements 
is < [Fe/H] >~ —1.5 dex, but that of the corrected data 
is < [Fe/H] >c± —1.0 dex. Also metallicity versus redshift 
studies are affected. For instance, a linear regression of 
[Fe/H] versus z a b s for the sample of Table 3 yields slopes^ 
to ~ —0.32 and ~ —0.17 for the data with and with- 
out dust correction, respectively. The slope is steeper for 
the corrected data, as expected if dust depletion is more 
important at higher metallicities. The slope m ~ —0.3 
of the corrected metallicities is similar to that found by 
Savaglio et al. (2000) and by Vladilo et al. (2000). In the 
former work a large sample of absorbers was corrected for 
dust effects assuming that the intrinsic abundance ratios 
in DLAs are solar. In the work by Vladilo et al. (2000) the 
uncorrected [Zn/H] ratio was used as a dust-free indica- 
tor of metallicity. While the interpretation of any possible 
trend between metallicity and redshift is very complex, it 

3 The slope is obtained from the least-squares method and 

. . „ , y".(^ta,i-<^abs»([Fe/H] i -<[Fe/H]» 

is defined as m — ! =^ A ^ , where 

x;,( z abs,i-<^b B » 2 

< ^abs > and < [Fe/H] > are the mean values of the sample. 
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-1.8 -<D.e 

[Fe/H] 

Fig. 2. [Si/Fe] and [Fe/H] ratios in DLAs. Crosses: orig- 
inal measurements. Dotted line: linear regression through 
the original measurements. Filled circles: measurements 
corrected with the set of parameters Ell (Table 2). 
Dashed line: linear regression through the dust-corrected 
data. Open squares: [S/Zn] versus [Zn/H] data in DLAs. 
Continuous line: midmean vector of [Si/Fe] versus [Fe/H] 
in Galactic metal-poor stars (Ryan et al. 1996). 

is clear that studies of this type based on [Fe/H] data not 
corrected for dust should be considered with caution. 



3.2.2. [Si/Fe] abundance ratios 

The abundance ratio of a-capture over iron-peak elements 
is an indicator of chemical evolution. In Galactic stars of 
low metallicity the a/Fe ratio is enhanced relative to the 
solar value and decreases with increasing [Fe/H] (e.g. Ryan 
et al. 1996). This behaviour can be naturally explained 
by models of galactic evolution in terms of a "time delay" 
between the enrichment of a elements and that of Fe-peak 
elements (e.g. Matteucci 1991 and refs. therein). 

Observed trends of a/Fe ratios in DLAs, particularly 
the Si/Fe ratio, have been used to characterize the evolu- 
tionary status of DLA galaxies (Lu et al. 1996; Prochaska 
& Wolfe 1999, 2002; Pettini et al. 2000; Molaro et al. 
2000). The [Si/Fe] ratios taken at face value are signifi- 
cantly enhanced at low metallicity (~ +0.4 dex), and this 
fact was originally taken as an indication that DLAs may 
be undergoing a chemical evolution similar to that of a 
young Milky Way (Lu et al. 1996). Subsequent work on 
dust free a/Fe ratios (Centurion et al. 2000; Molaro et al. 
2000) and on dust-corrected Si/Fe ratios (Paper I) has se- 
riously questioned this point of view. The enlargement of 
the abundance data base gives now a convincing evidence 
that the [Si/Fe] enhancement is largely due to dust deple- 
tion, as shown for instance by the behaviour of the [Si/Fe] 
versus [Zn/Fe] ratios (Prochaska & Wolfe 2002). In any 



15 -0.25 
0.0 0.0.0 



2.0.0 
redshift 
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Fig. 3. [Si/Fe] ratios versus absorption redshift in DLAs. 
Crosses: original measurements. Dotted line: linear regres- 
sion through the original measurements. Filled circles: 
measurements corrected with the set of parameters Ell 
(Table 2). Dashed line: linear regression through the dust- 
corrected data. Open squares: [S/Zn] ratios versus absorp- 
tion redshift in DLAs. 



case, as we show in Fig. ||, both the original [Si/Fe] mea- 
surements and the corrected ones indicate a behaviour of 
DLAs different from that typical of metal-poor Galactic 
stars. In fact, most of the non-corrected data lie above the 
"midmean vector" representative of Galactic stars (contin- 
uous line; Ryan et al. 1996), while most of the corrected 
data lie below. This is true for all the different sets of input 
parameters considered in Table 2. The filled circles in Fig. 
H indicate the data corrected with the set of parameters 
Ell. Very similar results are obtained for the other sets of 
parameters with enhanced [Zn/Fe]j ratio (E00, E01, E10). 
The corrected data lie even more below the Galactic stars 
adopting the sets of parameters with solar [Zn/Fe]j ratio 
(S00, S01, S10, Sll). 

The amount of [Si/Fe] enhancement at the lowest DLA 
metallicities (< —1.5 dex) is a probe of the earliest stages 
of the chemical evolution. From the analysis of the [Si/Fe] 
versus [Si/H] data, Prochaska & Wolfe (2002) find evi- 
dence for a [Si/Fe] plateau ([Si/Fe] ~ +0.3 dex) at [Si/H] 
< —1.5 dex, that they attribute to the lack of dust de- 
pletion at low metallicities. The present analysis indicates 
that 2 out of 6 DLAs with metallicities < — 1.5 dex have 
non-negligible [Si/Fe] depletion (Table 3). This indicates 
that some systems may be depleted even in the Prochaska 
& Wolfe's [Si/Fe] plateau. 

In addition to the degree of enhancement, the possi- 
ble presence of a slope in the [Si/Fe] data gives impor- 
tant information on the chemical evolution. Independent 
of the specific prescriptions of the models, the time delay 
between the injection of a elements and that of Fe-peak el- 
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ements implies that the a/Fe ratio must decrease with in- 
creasing [Fe/H]. The [Si/Fe] measurements in DLAs taken 
at face value do not show such a trend. The untilted dotted 
line in Fig. ^ represents the linear regression through the 
non-corrected data. On the other hand, the corrected ra- 
tios do show evidence for a decrease of [Si/Fe] with [Fe/H]. 
As an example, the tilted dashed line in Fig. || represents 
the linear regression through the data corrected with the 
Ell parameters. The slope resembles that of metal-poor 
stars, even if the ratios are lower on the average. 

Once corrected for dust depletion, the [Si/Fe] ratios 
show a tendency to increase with redshift which is not 
seen in the original measurements (Fig. ||). The typical 
value is [Si/Fe] w 0.0 dex at z — 1 and increases up to 
[Si/Fe] w +0.25 dex at z ~ 3.5, with a large dispersion at 
all redshifts. This possible trend of decreasing [Si/Fe] with 
cosmic time is consistent with the general expectations of 
the time delay models. Since DLAs probably represent an 
heterogeneous collection of different types of galaxies, the 
interpretation of this result is quite complex and will be 
considered in a separate work. In any case, the fact that 
only the dust-corrected ratios show the expected trends is 
rather encouraging. 

An important support to the validity of the correc- 
tion procedure comes from the comparison with [S/Zn] 
data. Since S and Zn are almost unaffected by deple- 
tion, the S/Zn ratio can be considered a sort of dust free 
a/Fe ratio (Centurion et al. 2000). If the method pre- 
sented here works properly, we expect that the corrected 
[Si/Fe] ratios and the non-corrected [S/Zn] ratios show 
a similar behaviour^] The limited sample of [S/Zn] data 
(Lu et al. 1996; Kulkarni et al. 1996; Molaro et al. 1998; 
Prochaska & Wolfe 1999; Centurion et al. 2000) shown as 
open squares in Figs. || and ||, confirms this expectation. 
The non-corrected [Si/Fe] ratios are generally inconsistent 
with the [S/Zn] data. 



4. Discussion and summary 

Chemical abundances of DLAs probe the effects of nucle- 
osynthesis in high H I density regions of the universe seen 
at different cosmic epochs throughout a considerable frac- 
tion of the Hubble time. However, disentangling the effects 
of nucleosynthesis from those of dust depletion is a rather 
difficult task. In this work a new method has been pre- 
sented for quantifying dust effects, with the aim of recov- 
ering the intrinsic abundances starting from the observed 
ones. There are several important differences with respect 
to the original method presented in Paper I and previous 
work on the same subject. 

The first difference concerns the definition of the dust- 
to-metals ratio. In Paper I this was defined in terms of 
extinction per unit level of metallicity, with the aim of 

4 Apart from little differences that may be present if S does 
not track perfectly Si, and Zn does not follow exactly Fe. 



estimating indirectly the extinction of DLAs from the el- 
emental depletions measured in each system. The broad 
agreement of such indirect estimates with the direct, but 
rather uncertain, estimates obtained from QSO reddening 
studies (Pei et al. 1991) is a positive test of consistency 
of the results obtained in Paper I. However, in order to 
compare these two types of estimates it is necessary to as- 
sume that the extinction is proportional to the amount of 
dust present in the system. This assumption is risky, since 
the extinction is also determined by the geometry and size 
distribution of the dust grains. The present work focuses 
only on the abundances in the dust and in the medium 
(gas plus dust). In order to avoid any assumption on the 
extinction of the dust, we have defined the dust-to-metals 
ratio as the fraction of atoms in dust form of an element 
Y chosen as a reference, i.e. r = Ny, dust/A/y, medium- 

The second difference is the way that elemental deple- 
tions are scaled up or down as a function of the dust-to- 
metals ratio. In the work by Kulkarni et al. (1997) and 
in Paper I the dust-to-metals ratio was allowed to change, 
but the abundances in the dust were kept constant. In this 
case, the fractions in dust of all elements are scaled up or 
down together. Therefore, all fractions in dust are propor- 
tional to that of the reference element, i.e. fx oc r. In the 
present work the fraction in dust of various elements have 
different dependence on r, according to the general law 
fx oc r 1+T,x . This law has been derived and calibrated in 
Paper II and is able to reproduce all the types of deple- 
tion patterns observed in the Milky Way interstellar gas 
(Savage & Sembach 1996) with a single set of rjx param- 
eters. 

The third difference is that, for the first time, elemental 
depletions are allowed to vary also as a function of the 
overall abundances of the medium. As shown in Paper II, 
the dependence of the fraction in dust on the intrinsic 
abundances must be of the type fx oc 10^ ex_1 ^W. The 
introduction of the parameter ex allows us to consider 
different behaviours of the dust composition. When ex = 
0, the composition of the dust does not vary with the 
composition of the medium. This extreme case could be 
valid, for instance, if the number ratio of two elements in 
the dust is uniquely determined by the capability of these 
two elements to have chemical bonds. In general, however, 
we may expect that the abundance ratio X/Y in the dust 
will track the abundance ratio X/Y in the medium, in 
which case ex = 1- By considering both cases ex = and 
ex = 1 we can take into account very different behaviours 
of the dust composition with respect to changes of the 
overall abundances, without entering into the details of the 
complex physical processes of dust formation, accretion 
and destruction. 

As in Paper I, the present method can be applied to in- 
dividual DLA systems. If the quality of the observational 
data is good enough, the method can be used for study- 
ing the abundances of distinct velocity components inside 
a given DLAs. The main steps of the procedure for re- 
covering intrinsic abundances can be summarized as fol- 
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Table 3. [Fe/H] and [Si/Fe] abundance ratios in DLAs: observed and dust-corrected values' 



Identifier 




Sabs 




[Zn/H] oba 




[Fc/H] oba 




[Si/Fe] oba 


Refs" 




[Fe/H] 




[Fc/H] 




[Si/F 


'Co* 




[Si/Fe] c c or 






































soo & 


Sll 


E00 & 


Ell 


SOO & 


Sll 


E00 & Ell 


0000-263 


3 


.3901 


-2 


.05 


± 





.09 


-2 


.10 


± 





09 





26 


± 


0.04 


Mol&OO 


-2 


.05 


± 


0. 


10 


-2.15 ± 
-2.15 ± 






09 
.09 


0.15 
0.15 


± 

± 


0.07 
0.07 






.25 ± 0.05 
.25 ± 0.05 


0149+33 


2 


.1400 


-1 


.65 


± 





.14 


-1 


.78 


± 





.10 





30 


± 


0.05 


P&W99 


-1 


.65 


± 


0. 


15 


-1.75 ± 
-1.75 ± 






12 
12 


0.18 
0.18 


± 

± 


0.12 
0.12 






.28 ± 0.08 
.28 ± 0.08 


0201+365 


2 


.4620 


-0 


.27 


± 





.00 


-0 


.88 


± 





04 





48 


± 


0.01 


P&W96 


-0 


.09 


± 


0. 


10 


-0.30 ± 
-0.27 ± 






08 
09 


-0.07 
-0.11 


± 
± 


0.01 
0.21 


-0 
-0 


.01 ± 0.03 
.03 ± 0.05 


0302-223 


1 


.0095 


-0 


.56 


± 





.12 


-1 


.20 


± 





12 





47 


± 


0.06 


Pet&OO 


-0 


.36 


± 


0. 


17 


-0.58 ± 
-0.55 ± 






.15 
.16 


-0.08 
-0.14 


± 
± 


0.03 
0.35 


-0 
-0 


.04 ± 0.06 
.07 ± 0.12 


0347-383 


3 


.0250 


-1 


.66 


± 





.11 


-1 


.81 


± 





.12 





.47 


± 


0.07 


Pro&Ol 


-1 


.06 


± 


0. 


15 


-1.76 ± 
-1.76 ± 






.14 
.14 


0.32 
0.32 


± 
± 


0.11 
0.11 






.42 ± 0.10 
.42 ± 0.10 


0454+039 





.8597 


-1 


.01 


± 





.11 


-1 


03 


± 





07 





.24 


± 


0.10 


Pet&OO 


-1 


.01 


± 


0. 


on 


-1.03 ± 
-1.03 ± 






07 
07 


0.22 
0.22 


± 

± 


0.11 
0.11 






.24 ± 0.10 
.24 ± 0.10 


0515-4414 


1 


.1510 


-0 


.99 


± 





.11 


-1 


.72 


± 





22 





46 


± 


0.27 


dlV&OO 


-0 


.70 


± 


0. 


44 


-0.97 ± 
-0.90 ± 






.37 
44 


-0.08 
-0.20 


± 

± 


0.09 
2.52 


-0 
-0 


.08 ± 0.17 
.15 ± 1.08 


0528-2505 


2 


.8110 


-0 


.76 


± 





.12 


-1 


.20 


± 


0. 


.15 





51 


± 


0.12 


Lu&96 


-0 


.67 


± 


0. 


20 


-0.83 ± 
-0.82 ± 






22 
23 


0.01 
0.01 


± 

± 


0.14 
0.15 






.11 ± 0.17 
.11 ± 0.13 


0841+129 


2 


.3745 


-1 


.49 


± 





.10 


-1 


.63 


± 





.17 





.37 


± 


0.15 


P&W99 


-1 


.49 


± 


0. 


23 


-1.59 ± 
-1.59 ± 






20 
20 


0.23 
0.23 


± 
± 


0.22 
0.22 






.33 ± 0.18 
.33 ± 0.18 


1104-1805 


1 


.6616 


-1 


.02 


± 





.01 


-1 


.59 


± 


0. 


02 





.64 


± 


0.03 


Lop&99 


-0 


.88 


± 


0. 


05 


-1.07 ± 
-1.05 ± 






04 
04 


0.06 
0.08 


± 
± 


0.03 
0.06 






.16 ± 0.03 
.17 ± 0.03 


1117-1329 


3 


.3511 


-1 


.24 


± 





.10 


-1 


53 


± 





09 





.26 


± 


0.06 


Per&02 


-1 


.23 


± 


0. 


13 


-1.34 ± 
-1.34 ± 






12 
12 


-0.03 
-0.03 


± 

± 


0.10 
0.10 






.07 ± 0.10 
.07 ± 0.10 


1210+1731 


1 


.8918 


-0 


.88 


± 





.10 


-1 


.14 


± 





.12 





.27 


± 


0.07 


Pro&Ol 


-0 


.87 


± 


0. 


14 


-0.98 ± 
-0.98 ± 






.14 
.14 


0.01 
0.01 


± 
± 


0.10 
0.09 






.11 ± 0.10 
.11 ± 0.10 


1215+33 


1 


.9990 


-1 


.27 


± 





.08 


-1 


.71 


± 





09 





24 


± 


0.06 


P&W99 


-1 


.21 


± 


0. 


14 


-1.35 ± 
-1.35 ± 






.12 
12 


-0.18 
-0.20 


± 

± 


0.06 
0.06 


-0 
-0 


.10 ± 0.08 
.10 ± 0.08 


1223+178 


2 


.4661 


-1 


.60 


± 





.10 


-1 


.85 


± 





.10 





.27 


± 


0.02 


Pro&Ol 


-1 


.59 


± 


0. 


11 


-1.70 ± 
-1.70 ± 






11 
11 


0.02 
0.02 


± 
± 


0.04 
0.04 






.12 ± 0.04 
.12 ± 0.04 


1247+267 


1 


.2232 


-1 


.04 


± 





.26 


-1 


.41 


± 


0. 


.13 





28 


± 


0.16 


Pet&99 


-1 


.01 


± 


0. 


37 


-1.13 ± 
-1.13 ± 






32 
33 


-0.09 
-0.09 


± 

± 


0.27 
0.23 






.01 ± 0.31 
.01 ± 0.30 


1331+170 


1 


.7764 


-1 


.28 


± 





.05 


-2 


.07 


± 





04 





63 


± 


0.00 


P&W99 


-0 


.94 


± 


0. 


07 


-1.24 ± 
-1.14 ± 






06 
07 


-0.02 
-0.06 


± 
± 


0.00 
0.45 



-0 


.01 ± 0.01 
.01 ± 0.24 


1351+318 


1 


.1491 


-0 


.36 


± 





.10 


-1 


.00 


± 


0. 


13 





45 


± 


0.16 


Pet&99 


-0 


.16 


± 


0. 


33 


-0.38 ± 
-0.35 ± 






26 

.30 


-0.09 
-0.16 


± 
± 


0.08 
0.91 


-0 
-0 


.05 ± 0.15 
.09 ± 0.30 


1354+258 


1 


.4200 


-1 


.60 


± 





.14 


-2 


.02 


± 





11 





29 


± 


0.16 


Pet&99 


-1 


.55 


± 


0. 


25 


-1.69 ± 
-1.68 ± 






21 

22 


-0.12 
-0.13 


± 
± 


0.18 
0.17 


-0 
-0 


.03 ± 0.21 
.03 ± 0.21 


2206-199A 


1 


.9200 


-0 


.39 


± 





.10 


-0 


.84 


± 





.10 





.44 


± 


0.02 


P&W97 


-0 


.32 


± 


0. 


11 


-0.47 ± 
-0.46 ± 






.10 
11 


-0.02 
-0.02 


± 
± 


0.02 
0.02 






.08 ± 0.03 
.08 ± 0.02 


2230+025 


1 


.8642 


-0 


.70 


± 





.09 


-1 


.18 


± 


0. 


09 





43 


± 


0.02 


P&W99 


-0 


.63 


± 


0. 


10 


-0.78 ± 
-0.77 ± 






09 
.10 


-0.04 
-0.05 


± 
± 


0.03 
0.02 






.05 ± 0.04 
.05 ± 0.03 


2231-0015 


2 


.0662 


-0 


.86 


± 





.10 


-1 


.41 


± 


0. 


.12 





.55 


± 


0.07 


P&W99 


-0 


.74 


± 


0. 


17 


-0.92 ± 
-0.90 ± 






.15 
.16 


0.00 
0.00 


± 
± 


0.07 
0.15 






.10 ± 0.09 
.10 ± 0.07 


2359-0216 


2 


.0950 


-0 


.75 


± 





.10 


-1 


.67 


± 





.10 





90 


± 


0.03 


P&W99 


-0 


.31 


± 


0. 


12 


-0.64 ± 
-0.49 ± 






12 

.13 


0.05 
0.17 


± 
± 


0.01 
1.00 






.12 ± 0.03 
.20 ± 0.68 



a When necessary, the column densities of the literature have been revised using the oscillator strengths given by Wclty et al. (1999). 
b References for the Hi, Sill, Fen, and Zn II column densities, with the following exceptions: 

Hi in Q0000-26 and Q2231-0015 (Lu et al. 1996); Foil in Q0841 + 129 (Centurion et al. 2000); Zn II in Q0347-383 (Levshakov et al. 2002). 

c For each system, the values in the first row arc relative to the sets of parameters SOO or E00; the values in the second row to the sets Sll or Ell; 

the corrected [Fc/H] ratios are equal in the cases SOO and Sll; sec Table 2. 
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dlV&OO: de la Varga ct al. 2000; L6p&99: Lopez ct al. 1999; Lu&96: Lu et al. 1996; Mol&OO: Molaro et al. 2000; Per&02: Peroux et al. 2002; 
Pet&97: Pettini et al. 1997; Pet&99: Pettini et al. 1999; Pet&OO: Pcttini ct al. 2000; P&W96: Prochaska & Wolfe 1996; P&W97: Prochaska & 
Wolfe 1997; P&W99: Prochaska & Wolfe 1999; P&W00: Prochaska & Wolfe 2000; Pro&Ol: Prochaska et al. 2001. 



lows. First, the dust-to-metals ratio is derived with Eq. 
(||). The most natural choice for the [X/Y] ratio to be 
used in that equation is the [Zn/Fe] ratio. In practice, the 
dust-to-metals ratio is derived from the observed [Zn/Fc] 
value and an educated guess of the intrinsic [Zn/Fe] value. 
Since the equations of the new method take into account 
the dependence of the abundances in the dust on the abun- 
dances in the medium, possible changes of the intrinsic 
[Zn/Fe] ratio in DLAs can now be treated self-consistently. 
Therefore the unknown value of the intrinsic [Zn/Fe] ratio 
is now considered as an input parameter which, at vari- 
ance with Paper I, may also assume values different from 
zero. At this point the abundances of elements other than 
Zn can be recovered from the observed ones by means of 
Eqs. (||) and (||). The dust-corrected abundances can be 
obtained for different sets of input parameters. In this way, 



we can test the stability of the results for different values 
of the £x parameter and of the intrinsic [Zn/Fe] ratio. 

The new method has been applied to correct the few 
existing interstellar abundances measured with accuracy 
in a nearby galaxy, namely the SMC. The results of this 
exercise support the validity of the procedure. In fact, for 
all the abundance ratios considered (Fe/H, Mg/Fe, Si/Fe, 
S/Fe, Cr/Fe, Mn/Fe, and Ni/Fe) the corrected interstel- 
lar abundances are in agreement with the intrinsic SMC 
abundances (known from stellar data) . The method yields 
consistent results for the different lines of sight considered. 
The importance of the dust correction procedure is illus- 
trated by the fact that the non-corrected interstellar data 
would indicate an SMC metallicity -2.0 < [Fc/H] < -1.2, 
a strong enhancement of the [a/Fe] ratios, and a high de- 
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grcc of inhomogeneity of the SMC. All these indications 
are contradicted by stellar data. 

The successful results obtained for a low-metallicity 
galaxy such as the SMC are encouraging since the fi- 
nal goal of the method is recovering the abundances of 
DLA systems, also characterized by a low metallicity level. 
There are other reasons why we feel justified in applying 
the dust correction method to DLAs systems. In fact, the 
adopted scaling law (Eq. 0) can successfully model deple- 
tion patterns originated in a variety of interstellar envi- 
ronments, with different physical conditions (e.g. cold vs. 
warm), locations (e.g. Galactic halo vs. disk; SMC), and 
extinction properties (e.g., strong vs. weak 217.5 nm emis- 
sion bump). 

The results obtained by applying the present method 
to the sample of DLAs with available Hi, Sin, Fen, and 
Znll measurements can be summarized as follows. 

The mean metallicity of the sample is < [Fc/H] >~ 
— 1.0 dex for the dust-corrected data, significantly higher 
than that of the original data, < [Fe/H] >~ —1.5 dex. 
The difference between observed and corrected metallici- 
ties can be as high as 1 dex in some cases. Also the [Fe/H] 
versus redshift relation in DLAs is affected by dust effects, 
the slope being ~ —0.32 and ~ —0.17 for the data with 
and without dust correction, respectively. 

Also the Si/Fe ratio, the most frequently measured 
a/Fe ratio in DLAs, is significantly affected by depletion. 
Most of the dust-corrected [Si/Fe] ratios lie below the me- 
dian values representative of Galactic stars. The [Si/Fe] 
enhancement at low metallicities ([Fe/H] ~ —1.5/ — 2.0 
dex) is typically [Si/Fe] ~ +0.2 dex. This value is some- 
what lower than the [Si/Fe] ~ +0.3 dex plateau of the 
original [Si/Fe] data found by Prochaska & Wolfe (2002) 
at the same metallicities, suggesting that some depletion 
is present also in that plateau. 

Possible trends of [Si/Fe] versus [Fe/H] and versus red- 
shift are also affected by depletion effects. The corrected 
[Si/Fe] data decrease with increasing [Fe/H] and increase 
with z a b s . Both trends are qualitatively consistent with 
the expectations of the time delay between the injection 
of Q-capture elements and that of Fe-peak elements, pre- 
dicted by models of chemical evolution. The uncorrected 
[Si/Fe] ratios do not show such trends. 

An important test of the dust correction method is 
provided by the comparison between Si/Fe and S/Zn data. 
In fact, the S/Zn ratio is an indicator of the a/Fe ratio 
which does not require dust correction to be applied since 
both S and Zn are essentially undepleted in the interstellar 
gas. The [Si/Fe] ratios show a behaviour similar to that of 
the [S/Zn] ratios only if the corrections are applied. This 
result supports the validity of the method. 

A detailed interpretation of the intrinsic DLA abun- 
dances is quite complex since DLAs represent a heteroge- 
neous, and probably biased, collection of interstellar re- 
gions associated with different types of galaxies. A de- 
tailed comparison of the dust-corrected abundances with 
the predictions of chemical evolution models of different 



types of galaxies will be presented in a separate paper 
(Calura et al. 2002). 

Acknowledgements. I wish to thank Irina Agafonova for pro- 
viding routines for the solution of Eq. (^|). The data for the 
DLA system towards BR 1117-1329 have been kindly provided 
by Celine Peroux and collaborators in advance of publication. 
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